Plasmon bands in metallic nanostructures 
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The photonic band structure of a three-dimensional lattice of metal spheres is calculated using 
an embedding technique, in the frequency range of the Mie plasmons. For a small filling factor of 
the spheres, Maxwell-Garnett theory gives an almost exact description of the dipole modes, and 
the multipole modes are fairly dispersionless. For a larger filling factor, crystal field effects modify 
the multipole frequencies, which show dispersion. These multipole bands are enclosed between the 
dipole modes. For touching spheres, there is a wide continuum of plasmon modes between zero 
frequency and the bulk metal plasmon frequency, which yield strong absorption of incident light. 
These plasmon modes are responsible for the blackness of colloidal silver. 

PACS numbers: 42.25.Bs,42.70.Qs,73.22.Lp 



An isolated free-electron metallic sphere shows Mie 
plasmons at frequencies 
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where uj p is the bulk plasmon frequency, and / is the 
quantum number identifying the multipole. This assumes 
the Drude dielectric function, 
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where r is the scattering time. In a lattice, crystal 
field effects shift the Mie plasmon frequencies, an effect 
which has been studied both experimentally and theo- 
retically 0. These theoretical studies have generally 
been restricted to the dipole plasmons and their crys- 
tal field, but recently a number of methods have been 
developed for a full solution of Maxwell's equations in 
periodic structures 0, 0- El 0- ■ There have been pho- 
tonic band structure calculations of lattices of metallic 
rods 0, cylinders @, and spheres HE OH El- The 
plasmon modes and their interaction with light have been 
studied by Yannopapas et al. 0], in the low filling frac- 
tion regime. 

In this paper, we use the new embedding method to 
solve Maxwell's equations and consider the evolu- 
tion of plasmon bands in a lattice of nanoscale metallic 
spheres as the sphere size increases from very small, prac- 
tically the isolated sphere limit, to touching spheres. Our 
band structures and densities of states show clearly the 
interaction of light with such systems, and the way that 
the plasmons are affected by the crystal field. We find 
that the presence of structured metal introduces a contin- 
uum of plasmon modes, which yield strong absorption of 
incident light and are responsible for the unique optical 
properties of colloidal metals. 



In the photonic embedding method, the dielectric ob- 
jects - here metallic spheres - are replaced by an embed- 
ding "potential" over their surface. This is added to the 
wave equation for the uniform vacuum region between 
the spheres, and the equation is only solved explicitly in 
this region. The embedding potential, in reality a surface 
tensor operator, ensures that the surface parallel compo- 
nents of the E and H fields match across the boundary 
between the vacuum and the metal spheres. The electro- 
magnetic field is expanded in terms of any suitable basis 
set (we use vector plane waves) in the region between 
the spheres, and because this field has no discontinuities 
the expansion converges very well. In the case of non- 
metallic dielectric spheres, the convergence is at least an 
order of magnitude better than a direct plane wave ex- 
pansion through the whole of space (the usual method 
for finding photonic band structures) 13] ; in the difficult 
case of metallic spheres, the embedding method contin- 
ues to converge well. Of course there is a price to pay, 
and this is the evaluation of the embedding tensor which 
replaces the metallic spheres. This involves finding the 
exact solution of Maxwell's equations inside the spheres 
for each multipole / at frequency u>, but we have derived 
straightforward expressions that may readily be coded. 

The multipoles of the embedding tensor are cut off at 
a maximum value, Z m ax, imposing a limit to the plasmon 
modes on each sphere. Without this, the density of states 
increases without limit as lo approaches the planar sur- 
face plasmon frequency ui p / \[2 (the limit of large I in Eq. 
(1)). It is also necessary for convergence of the plane 
wave expansion. There is a physical plasmon cut-off in 
Z max due to Landau damping, though this is larger than 
the values we use. Controlling Z max also enables us to 
study the evolution of the plasmon bands in more de- 
tail. Unlike in scattering theory where Z max is similarly 
imposed, all higher values of I are included in our embed- 
ding, but with a different boundary condition, that the 
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FIG. 1: nk, spectral density for fee lattice of metal spheres, 
f = 2, uj p = 0.1, Z max =3, at k = (0.03,0,0) (solid line) and 
(0.1,0,0) (dashed line). 



FIG. 2: Band structure for fee lattice of metal spheres, f = 
0.5, cJ p = 0.1, in FX direction (solid lines and shaded area). 
The dashed line represents the free-space dispersion relation 
uj = k. Filled circles are results from Maxwell-Garnett theory. 



surface parallel components of H vanish on the surface 
of the sphere. 

The system we shall discuss here is a face-centred cubic 
(fee) lattice of metallic spheres in vacuum, the spheres 
being described by the Drude dielectric function of Eq. 
(2). The scaling property of Maxwell's equations makes 
it natural to use the dimensionless reduced frequency u> — 
aja/2nc and reduced wave- vector k = ka/27r, where a is 
the conventional lattice constant. Throughout this paper 
we use the reduced value of lu p — 0.1; taking uj p = 15 eV 
for Al, this corresponds to a lattice constant of a = 83 
A. The reduced lifetime f is taken to be 1000, a value 
which enables us to separate spectral features. We note 
that quantum size effects are likely to be important in 
the small spheres we study here - this will be treated in 
later work. 

In the embedding method, with a frequency-dependent 
embedding tensor, it is most straightforward to use a 
plane- wave expansion for the Green function between the 
spheres, and from this we calculate the spectral density 
n^Lu) (proportional to the energy density) integrated 
through the vacuum region. Figure ^ shows a typical 
spectral density, for an fee lattice of metal spheres with 
reduced radius f = 2irr/a = 2, corresponding to an ac- 
tual radius of r = 26. 3A and a filling fraction of 54%, at 
k = (0.03, 0, 0) and (0.1, 0, 0) (the X-point corresponds to 
k = (1, 0, 0)). In this calculation / max =3, and convergence 
has been achieved with 18 transverse and 181 longitudi- 
nal plane waves. The number of transverse waves stays 
very low over the frequency range shown, but the num- 
ber of longitudinal waves increases as the radius of the 
sphere decreases and / max increases. The large number of 
longitudinal waves required is characteristic of plasmon 
systems. The almost constant features of Fig. ^ between 
Co = 0.055 and 0.085 are the multipole plasmons with 
a non-degenerate dipole plasmon at the top of this fre- 



quency range. The dispersing peaks below and above this 
range are the light line coupled to a doubly degenerate 
dipole mode. 

From the peaks in the spectral density we plot the 
photonic band structure. We begin by considering an fee 
lattice of relatively small spheres with f = 0.5, an actual 
radius of 6.6 A and a filling fraction of 0.8%. The band 
structure for l max = 1 (solid lines in Fig. [2J) shows flat 
dipole modes at the Mie frequency of w p /V3 (solid hor- 
izontal line) and the doubly degenerate light line (solid 
lines below and above the flat dipole modes). At the 
crossing point the light lines mix with two of the dipole 
modes, but with small spheres the interaction and the 
splitting are small. This calculation required 2 trans- 
verse and 1591 longitudinal plane waves for convergence. 
Increasing Z max does not change the light line and flat 
dipole modes, but gives the flat multipole bands at the 
frequencies given by Eq. (1) filling the shaded area of 
Fig. El 

At low filling fractions, the interaction between the 
spheres can be considered within the Maxwell-Garnett 
approximation [l^j. For the light mode that we see in 
Fig. |21 we can define a frequency-dependent effective di- 
electric function 

e eS (cj) = ~k 2 /Q 2 , (3) 
and the Maxwell-Garnett approximation then gives (l5| 
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FIG. 3: Band structure for fee lattice of metal spheres, f — 
1.2, ujp = 0.1, i max — ► ex), in FX direction (solid lines and 
shaded area). Circles are results obtained with Z max = 1. 
Maxwell- Garnett calculations, not represented here, nearly 
coincide with l max — 1 results. 



FIG. 4: Band structure for fee lattice of metal spheres, 
f — 2.0, Up — 0.1, Z max = 2, in TX direction (solid lines). 
Filled circles are results from Maxwell-Garnett theory. 



Here / is the filling fraction, m = (1 — /)/3, and n — 
(1 + 2/)/3. The filled circles in Fig. H show the results 
for the dispersion calculated from Eqs. and (0J. There 
is, not surprisingly, precise agreement with our calculated 
results for the light dipole mode. 

The concept of effective dielectric function is particu- 
larly useful for understanding optical absorption, which 
is dictated by the poles of e c g. Conversely, the energy 
loss of fast charged particles is dictated by the poles of 
e^. An inspection of Eqs. (0J and (JSJ shows that in the 
isolated sphere limit (/ — > 0) both optical absorption and 
energy loss occur at the Mie dipole frequency of ui p /y/3. 
At finite filling fraction, the Maxwell-Garnett approxi- 
mation predicts optical absorption at y/mu) p and energy 



loss at y/nui p . These are, respectively, the limit of the 
lower branch of the light line and the starting frequency 
of the upper branch. 

At a sphere radius of f = 1.2, filling fraction 12%, 
we find that the plasmon frequencies deviate from the 
isolated sphere results. The circles of Fig. 3 show the 
band structure with Z max = 1. We see that the interaction 
with the light line produces a greater splitting of the 
dipole modes than in Fig. [21 which is indeed accurately 
reproduced by Maxwell-Garnett theory. Increasing Z max 
(solid lines of Fig. EJ) barely affects the lower and upper 
branch of the light line, only introducing virtually flat 
multipole modes at the Mie frequencies of Eq. which 
fill the shaded area of Fig. [3J 

A calculation of the band structure of a simple cubic 
lattice of Drude spheres with Q p — 0.06 and filling frac- 
tion of 12% was carried out by Pendry [9j, by using a 
transfer-matrix scheme. We have calculated the band 
structure for those parameters, and have found that it 
has almost exactly the same form as the band structure of 




FIG. 5: ilk, spectral density for fee lattice of metal spheres, 
f = 2, uj p = 0.1, at k = (0.5,0,0). Solid line, Z max = 2; long- 
dashed line, i max = 6; short-dashed line, i max = 8; dotted line, 
' max = 10; dashed-dotted line, Z max = 12. 



Fig. [3J This disagrees, however, with the band structure 
reported in Ref. 9], which shows flat low-energy bands 
absent in the band structure of Fig. [3J These probably 
originate from spurious modes associated with the edges 
and corners that are present in the transfer-matrix dis- 
cretization procedure. 

The plasmon band structure is much more interesting 
at f = 2, filling fraction 54%, with greater interaction 
between all the modes. Figure 0] exhibits the band struc- 
ture for Z max = 2, in the FX direction. We have also 
looked at this band structure in the TW direction, and 
have found that along this direction all the states are 
non-degenerate, with 8 plasmon bands (there is no / = 
plasmon). Figure 0] shows that the Maxwell-Garnett dis- 
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FIG. 6: Band structure for fee lattice of metal spheres, 
f = 2.0, uj p = 0.1, Z max — » oo, in FX direction (solid lines and 
shaded area). The circles show the dipole bands of Fig. 2] 
obtained with / max = 2. 




FIG. 7: Band structure for fee lattice of touching metal 
spheres, f — 2.2, lo v = 0.1, 2 max = 12, in TX direction. Filled 
circles are results from Maxwell-Garnett theory. 

persion calculated from Eqs. and (0J starts to become 
somewhat inaccurate, as its assumptions begin to break 
down. We also note that unlike in the case of smaller 
filling fractions, where the high multipole plasmons cut 
across the upper branch of the light line and lie higher 
than the flat dipole plasmon modes, here the dipole plas- 
mons enclose the multipole modes. 

As ? max increases, there are so many plasmon bands 
that they become difficult to track individually. How- 
ever, we can see the general behaviour by studying the 
spectral density. Figure shows the spectral density at 
k = (0.5, 0, 0) for f — 2 and various values of Z max : 2, 6, 
8, 10, and 12. We see that there is very little difference 
in the overall structure for Z max > 6, the only difference 



being the large increase of the spectral density associated 
with high-/ multipole plasmons close to the planar sur- 
face plasmon limit at oJ p /v2 = 0.07, in the middle of the 
multipole plasmon region. We can then plot the band 
structure in the large Z max limit, Fig. (convergence has 
been achieved with Z max = 8), where the multipole bands 
lie in the shaded area enclosed by the dipole bands. The 
dipole bands themselves have been pushed out compared 
with i max = 2 (shown in Fig. by circles). 

It is clear from our results for f = 2, showing an in- 
crease in dipole band gap as / max increases, that part of 
the band gap results from the interaction of the dipole 
modes with the multipole modes. This interaction means 
that the multipole modes can be excited by both photons 
and by fast charged particles. This contrasts with the ap- 
parent absence of interaction in Fig.|3for f = 1.2. 

The behaviour of the multipoles is more extreme once 
the spheres touch 0, corresponding to f = 2.2, and a 
filling fraction of 74%. Figure shows the band struc- 
ture with Z max = 12. The multipole modes in the shaded 
area are still bounded by the dipole modes, but as Z max 
increases the multipole modes steadily broaden, pushing 
the lower dipole branch towards Co — and the upper 
branch towards uj p . The upper branch appears to have 
reached its upper limit by Z max = 12, just below oJ p ; the 
lower branch is still moving down, and it is not clear 
whether this reaches the zero frequency of the bulk metal. 
However, there is a very wide band of multipole modes 
nearly spanning these limits, with corresponding optical 
absorption. This is the reason for the striking blackness 
of colloidal silver. 
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